Anisotropic itinerant magnetism and spin fluctuations in BaFe2As2: A neutron 

scattering study 
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Neutron scattering measurements were performed to investigate magnetic excitations in a single- 
crystal sample of the ternary iron arsenide BaFe2As2, a parent compound of a recently discovered 
family of Fe-based superconductors. In the ordered state, we observe low energy spin-wave exci- 
tations with a gap energy A = 9.8(4) meV. The in-plane spin-wave velocity Vat and out-of-plane 
spin-wave velocity Vc measured at 12 meV are 280(150) and 57(7) meV A, respectively. At high 
energy, we observe anisotropic scattering centered at the antiferromagnetic wave vectors. This scat- 
tering indicates two-dimensional spin dynamics, which possibly exist inside the Stoner continuum. 
At Tjv = 136(1) K, the gap closes, and quasi-elastic scattering is observed above Tjv, indicative 
of short-range spin fluctuations. In the paramagnetic state, the scattering intensity along the L 
direction becomes "rodlike," characteristic of uncorrelated out-of-plane spins, attesting to the two- 
dimensionality of the system. 

PACS numbers: 75.30.Ds, 78.70.Nx, 74.72.-h,75.50.Ee 



I. INTRODUCTION 

A pairing mechanism mediated by magnetic fluctua- 
tions has been regarded as a leading candidate to re- 
solve the problem of high- Tc superconductivity.- The 
recent discovery of superconductivity in a family of 
iron arsenideSj^ whose non-superconducting parent com- 
pounds exhibit long-range antiferromagnetic order sim- 
ilar to the cuprates^ provides yet another promising 
group of materials for studying the interplay between 
magnetism and superconductivity. One property shared 
by most high-Tc superconductors is two-dimensional 
magnetism, and both cuprates and iron arsenides are 
comprised of magnetic layers, sandwiched between non- 
magnetic ions. However, while parent compounds of the 
cuprates are Mott insulators with large in-plane exchange 
interactions (J = 135 meV for La2Cu04),4 magnetism in 
the metallic iron arsenides most likely originates from 
itinerant electrons, and the antiferromagnetic order is 
a result of a spin-density-wave (SDW) instability due 
to Fermi-surface nestingi^ In order to study anisotropic 
magnetic interactions and itinerant magnetism in the 
iron arsenides, measurements of magnetic excitations on 
a single-crystal sample are extremely crucial, and can 
provide new insights into the understanding of high-Tc 
superconductivity. 

Shortly after the discovery of superconductivity in 
fluorine-doped LaFeAsO)^ Rotter et al.^ proposed the 
ternary iron arsenide BaFe2As2 as a "very promising 
candidate" for superconductivity based on its almost 
identical structure and electronic properties. Similar to 
LaFeAsO, BaFe2As2 is comprised of FeAs layers, and 
shows a transition to an antiferromagnetic ordered state 
and a structural transition from the tetragonal space 
group lA/mmm at high temperature to the orthorhombic 
space group Fmmm at low temperature; though, those 



transitions appear at different temperatures in the for- 
mer but at the same temperature in the latteri^ Sub- 
sequently, several groups reported superconductivity in 
potassium, cobalt, cesium, and nickel doped AFe2As2, 
where A—Ba., Sr, Ca, and Eu, with the highest Tc of 38 
K in Bao.6Ko.4Fe2As2 4^'iaaii2 

Powder neutron diffraction on BaFe2As2 reveals the 
long-range magnetic order below the Neel temperature 
Tn = 136(1) K, where Fe spins with an ordered moment 
of 0.87(3) fiB per Fe form a colinear spin structure along 
the crystallographic a axis, and antiferromagnetic (ferro- 
magnetic) arrangement along the a and c axes (6 axis)ii 
Previous inelastic neutron scattering studies on a pow- 
der sample of BaFe2 As2 show magnetic excitations possi- 
bly up to 170 meV, indicative of strong spin couphngs»i^ 
Neutron scattering studies on single-crystal samples of 
the related compounds SrFe2As2 (T^r = 220 K) and 
CaFe2As2 (T/v = 172 K) reveal low energy spin- wave 
excitations with spin gaps of < 6.5 and 6.9(2) meV, 
respectively] ^"^i Above the ordering temperature, the 
spin gaps close, and the spin-wave scattering is replaced 
by broad quasi-elastic scattering. ^5 

Here, we present the flrst inelastic neutron scattering 
study on a single-crystal sample of BaFe2As2. We ob- 
serve quasi-elastic scattering, which peaks in the vicinity 
of Tjv, indicative of short-range spin fluctuations. Above 
Tjv, the out-of-plane spins become uncorrelated and scat- 
tering intensity becomes "rodlike" along L, revealing the 
two-dimensionality of the system. Below Tjv, we observe 
low energy spin-wave excitations with a spin gap. The 
difference between in-plane and out-of-plane spin-wave 
velocities suggests strongly anisotropic magnetism. Most 
interestingly, at high energy, we observe anisotropic scat- 
tering centered at the antiferromagnetic wave vectors. 
This scattering possibly indicates magnetic excitations 
inside the Stoner continuum. 
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FIG. 1: (Color online) (a) Constant-energy scans at Two = 0.7 
meV and T = 136 K measured in the flat mode on Cl-1. 
Solid line shows a fit to a Gaussian convoluted with the in- 
strumental resolution function and the dotted line shows the 
instrumental resolution, (b) Temperature dependence of the 
magnetic Bragg intensity measured at (1,0,3) and of magnetic 
fluctuations measured at Q = (1, 0, 1) and fiuj = 0.7 meV us- 
ing the seven-blade horizontally focused analyzer. The line 
corresponds to a power- law form Af^^ « [Tn - Tf^ for sub- 
lattice magnetization (which is proportional to the square 
root of the Bragg intensity) with the exponent /? — 0.125 and 
Tn = 136(1) K. (c) The contour plot as a function energy 
and Q shows scattering intensity around (1,0,1) at T = 136 
K using the three-blade horizontally focused analyzer. 



II. EXPERIMENTAL DETAILS 

A single-crystal sample, which is composed of two 
coaligned crystals of total mass 0.4 g, was grown us- 
ing the Bridgman method with the FeAs flux. The 
detailed synthesis and sample characterizations are re- 
ported elsewhereJi Neutron-diffraction measurements 
on this sample yield the ordered moment of 0.91(0.21) /is 
per Fe, which is consistent with the previous powder- 
diffraction work.^ Inelastic neutron scattering measure- 
ments were performed on the triple-axis spectrometers 
HER (Cl-1) (cold neutrons) and GP-TAS (4G) (ther- 
mal neutrons), which are operated by the Institute for 
Solid State Physics, University of Tokyo. The sample 
was aligned in the hOl zone with the reciprocal-lattice pa- 
rameters a* = 1.124(3) and c* = 0.485(2) A^^. The 
final energy was fixed at 5 meV at Cl-1 and 14.7 meV at 
4G. Vertically focused pyrolytic graphite (PG) crystals 
were used to monochromate the incident neutron beam 
using the 002 reflection. Vertically focused (horizontally 
flat mode) and doubly focused PG crystals were used to 
analyze the scattered-neutron beam at Cl-1, while ver- 
tically focused PG crystals were utilized to analyze the 
scattered-neutron beam at 4G. At Cl-1, horizontal coUi- 
mations of guide— open— sample— 80' — 80' were employed 
for the flat mode, and guide— open— sample— radial— open 



for the horizontally focusing mode. At 4G, horizontal col- 
limations of 40' — 80'— sample— 80' — 80' were employed. 
PG filters were placed in the scattered beam at 4G, and 
cooled Be/oriented-PG crystals and room-temperature 
Be filters were placed in the incident and scattered beams 
at Cl-1, respectively, to reduce higher-order contamina- 
tion. The sample was cooled by a closed cycle ^He cryo- 
stat. 



III. SPIN FLUCTUATIONS NEAR Tn 

A constant-energy scan around (1,0,1) at huj = 0.7 
meV and T = 136 K performed on Cl-1 using the hori- 
zontally flat analyzer [Fig. [TJa)] shows spin fluctuations 
in the vicinity of T/y. The temperature dependence of 
the low energy scattering at Q = (1,0,1) and fiuj = 0.7 
meV using the seven-blade horizontally focused analyzer 
[Fig. [TJb)] shows a peak at T/v, which coincides with 
the onset of the magnetic Bragg intensity measured at 
(1,0,3). As a function of temperature, this scattering 
intensity rapidly decreases on the low-temperature side, 
resulting from opening of the spin gap (the detailed dis- 
cussion of the spin gap will be presented below), but 
has much weaker temperature dependence on the high- 
temperature side, indicative of critical scattering above 
Tjsf. Several constant-energy scans performed at 136 K 
using the three-blade horizontally focused analyzer con- 
stitute a contour plot [Fig. [ijc)], which shows a non- 
dispersive scattering column centered at the antiferro- 
magnetic Brillouin-zone center (1,0,1) that appears to 
extend up to high energy. The scattering column, which 
exists at least up to 2 meV (the highest energy of the Cl- 
1 measurements), was also observed around (1,0,3) at 136 
K (not shown). The constant-energy scans at hoj = 0.7 
meV measured at Cl-1 [Fig. [ija)] and at Tiuj = 12 meV 
measured at 4G (not shown) yield correlation lengths 
^ = 15(1) A and ^ = 18(2) A, respectively, corresponding 
to roughly 6 times the nearest-neighbor distance. 

Our neutron scattering measurements show that the 
scattering near T/v are dominated by short-range spin 
fluctuations. The presence of these short-range fluctu- 
ations might appear to contradict a report of the first- 
order magnetic transition observed in the ''^As nuclear- 
magnetic-resonance study.™ However, large magnetic in- 
teractions can give rise to the short-range spin fluctua- 
tions in the paramagnetic state, despite the first-order 
magnetic transition. In addition, two-dimensionality in 
this layered system can enhance spin fluctuations. The 
temperature dependence of the order parameter mea- 
sured with neutron scattering shows no sign of hystere- 
sis. Furthermore, the spin gap decreases continuously 
as a function of temperature (not shown), showing no 
sign of an abrupt change at T^r that is expected for a 
first-order transition. On the other hand, we do not ob- 
serve the divergence of correlation length at Tn, which 
would suggest a second-order transition. Therefore, the 
nature of the magnetic transition in BaFe2As2 cannot be 
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FIG. 2; (Color online) Constant-energy scans around (1,0,1) 
at hcu = 10 and 12 meV, and at T = 5 and 145 K. (a) -(d) 
show the scans along H and (e) — (h) along L. Solid lines in 
(b), (d), (f), and (h) show fits to the empirical dispersion re- 
lation described in the text convoluted with the instrumental 
resolution function. Solid lines represent the best fits with 
A = 9.8(4) meV, Vat = 280(150) meV A, and Vc = 57(7) 
meV A. The L scans at T = 145 K, (a) and (b), are fit with 
Gaussians convoluted with the resolution function. Dotted 
lines in (a) — (d) indicate the instrumental resolution, while 
the dotted lines in (e) and (g) show the background level es- 
timated from the low-temperature data. 



resolved by our neutron scattering measurements. How- 
ever, neutron scattering measurements of the order pa- 
rameter on a related compound a-FeTe show thermal 
hysteresis at the transition, supporting the first-order 
magnetic transitions^ 



IV. LOW ENERGY SPIN- WAVE EXCITATIONS 

In the ordered state, scattering from low energy spin- 
wave excitations was observed. Figures [SJ^b), 2(d), 2(f), 
and 2(g) show a series of constant-energy scans around 
(1,0,1) measured at 4G at T = 5 K. The observed peaks 
were fit with narrow (resolution-limited) Lorentzians 

L^{q,huj) = ^^^"^ inry^M^^-n^^y' convoluted with the 
instrumental resolution function taking into account the 



following empirical spin-wave dispersion relation: 

?i^q - ,jA^+vl,{ql+4) + (^cq^f, (1) 

where A is a gap energy in meV, and Vat and Vc are in- 
plane and out-of-plane spin-wave velocities in a unit of 
meV A, respectively, q — {qx , Qy , q^) is a wave vector in 
the first Brillouin zone away from the zone center Qa in 
a unit of A~^. It should be noted that this spin- wave 
model is used as a tool to analyze the data, and does not 
imply a localized spin model. A was obtained by fitting 
the 10 meV constant-energy scans, since these scans were 
performed in the proximity of the spin gap and, therefore, 
are the most sensitive to A. Vab and Vc were obtained by 
fitting the 12 meV constant-energy scans. Using the fit 
parameters from the previous round of fitting as inputs, 
we repeated the fitting several times to assure that all fit 
parameters are self-consistent, since the three parameters 
are not independent due to a finite instrumental resolu- 
tion. The obtained fit parameters are A = 9.8(4) meV, 
Vab = 280(150) meV A, and Vc = 57(7) meV A, where 
the error corresponds to three times the statistical error. 
The gap energy is consistent with a constant-Q scan mea- 
sured at (1,0,3) [Fig. [3]Ja)]. The overall prefactor A{huj) 
follows the l/fku scaling within the error bars, consistent 
with a report in Rcf. 15. The black solid lines in Fig. [2] 
show the best fits to the empirical dispersion relation 
[Eq. [1] using the aforementioned parameters. The line 
shape and apparent shift from the Brillouin-zone center 
are governed by the four-dimensional resolution function. 
The dotted lines indicate the instrumental resolution, as- 
suming a (5-function excitation. The spin-wave velocities 
are expected to increase at high energy away form the 
spm gap; Fig. Hi shows Vc = 80(10) meV A at ?iw = 24 
meV. Unfortunately, as will be discussed below, a scan 
along H at high energy does not allow accurate determi- 
nation of Vab- 

As a comparison, the spin-wave velocities for CaFe2As2 
are Vab = 420(70) meV A and Vc = 270(100) meV A,^^ 
and for SrFe2As2, Vab = 560 meV A and Vc = 280 
meV A (note a factor-of-2 difference from Ref. [141 ). Com- 
pared to CaFe2As2 and SrFe2As2, the much smaller Vc 
suggests the significantly weaker interlayer coupling, and 
hence lower Tjv in BaFe2As2. Furthermore, the smaller 
ratio of Vc to Vab indicates the more anisotropic mag- 
netic interactions in this system. As a further compar- 
ison, the spin gap in BaFe2As2 is larger than those in 
CaFe2As2 and SrFe2As2. The origin of the spin gaps in 
this family of FeAs compounds is still under much de- 
bate. In other itinerant antiferromagnetic systems, for 
example, in a family of 7-Mn alloys, where spin gaps 
of 7-10 meV were observedji2i2£i2i Fishman and Liu3^ 
have successfully derived the temperature dependence of 
the spin gaps in agreement with the experiments using 
a two-band model with magnetoelastic interactions. In 
the iron arsenide systems, the closeness of the structural 
and magnetic phase transitions, which has also been ob- 
served in the 7-Mn alloys, provides tantalizing evidence 
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for strong magnetoelastic interactions that could explain 
the presence of the spin gaps. In addition, dipole inter- 
actions and spin-orbit couplings could also contribute to 
the presence of the spin gap in these systems. 

At T = 145 K above T/v, the spin- wave scattering is re- 
placed by broad quasi-elastic scattering that extends up 
to high energy. The H scans in Figs. Ufa) and 2(c) yield 
correlation lengths ^ = 13(2) A at ?ia; = 10 meV, and 
f = 15(2) A at = 12 meV, which are slightly smaller 
at this temperature compared to 136 K as the correla- 
tion lengths are expected to decrease as temperature in- 
creases. On the other hand, the L scans in Figs.[2je) and 
2(g) show rodlike scattering, indicative of uncorrelated 
interlayer spins. The two-dimensionality of this system is 
clearly observed in the paramagnetic state with large in- 
plane correlation lengths and uncorrelated out-of-plane 
spins. Similar to the cuprates, the relatively large cor- 
relation lengths in the paramagnetic state are a conse- 
quence of the large magnetic interactions and the two- 
dimensionality of the magnetismi^ Despite this similarity, 
however, the origin of magnetism in the iron arsenides is 
believed to be quite different from the cuprates. 



V. MAGNETIC EXCITATIONS AT HIGH 
ENERGY 

In order to further investigate the nature of magnetism 
in BaFe2As2, we measure inelastic scattering in an abso- 
lute unit of barn/meV at T = 5 K to estimate a fluc- 
tuating moment. A normalization function to inelastic 
scattering was obtained from acoustic phonons measured 
on copper. The normalized intensity can be described by 
the following equation: 

/(q, hi,) = (^) ' [f{Q)f s{nLo)6{huj - huj^), (2) 

where (^)^ = 72.65 x 10^^ barn//^|, /(Q) is a mag- 
netic form factor for Fe^^, and ?iWq is the spin- wave 
energy described by the dispersion relation in Eq. [1] 
Figure [31[a) shows the scattering intensity S{hu!q) = 

(^^Y S{hujq) in barn/Fe measured at Q = (1,0,3) and 
(1,0,5) after corrected for background, absorption, the 
instrumental resolution, and the magnetic form factor. 
The same scattering was also observed at the antiferro- 
magnetic wave vectors (3,0,1) and (1,0,7), whose differ- 
ent intensities can be accounted for by the magnetic form 
factor and polarization effect. To obtain the fluctuating 
moment Mf, we calculate 

Mj = (——^ J S{huj)S{hLj - ULUq) dq doj. 

Since the ordered moment is along the crystallographic 
a axis, nearly equivalent intensities at (3,0,1) and (1,0,3) 
suggest that Mf are dominated by transverse fluctua- 
tions. Figure ^h) shows a product of S{hujq) and the 
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FIG. 3: (Color online) (a) Scattering intensity in barn/meV 
at T = 5 K measured at (1,0,3) (red circles) and (1,0,5) (blue 
squares) after corrected for background, absorption, the in- 
strumental resolution, and the magnetic form factor. Dotted 
line shows the 1/hui dependence. Shaded areas are guides to 
the eyes, representing the spin-wave scattering (yellow) and 
the "extra" scattering (gray), (b) shows a product of S{hu!ci) 
and the density of states. Solid line is a guide to the eyes. 
Inset shows the density of states. 



density of states, assuming the dispersion relation in 
Eq. [1] The integral was calculated up to 42 meV. We 
obtain Mf = 0.8(2) fj.B, which already accounts for most 
of the ordered moment. However, the magnetic excita- 
tions were reported to extend up to much higher energy^— 
This discrepancy suggests magnetic scattering from other 
origin. 

In fact, constant-energy scans at huj = 24 meV show 
an "extra" scattering intensity, which is centered at 
the antiferromagnetic wave vectors and coexisting with 
the spin-wave scattering. The gray shaded areas in 
Figs. IDJa)— 4(d) show this scattering around (1,0,5) and 
(3,0,1), which is sharp along H, but broad along L. 
The yellow shaded areas indicate the spin-wave scatter- 
ing with A = 9.8(4) meV, Vab = 280 meV A, and Vc = 80 
meV A; note a larger Vc at this energy. Unfortunately, the 
presence of the extra scattering in the H scans hinders ac- 
curate determination of Vat, which presumably becomes 
larger than 280 meV A at this energy. Consequently, the 
shaded areas in Figs. IDJa) and 4(c) are rough estimates 
and only serve as guides to the eyes. Figures IH^b) and 
4(d) describe the scattering intensities along L, which 
clearly show the coexistence of the spin-wave scattering 
and the extra scattering. The difference of the spin- wave 
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FIG. 4: (Color online) Constant-energy scans at 24 meV 
around [(a)-(b)] (1,0,5) and [(c)-(d)] (3,0,1). Solid lines 
show the spin-wave scattering with A = 9.8(4) meV, Vab ~ 
280 meV A, and Vc = 80 meV A. Dotted line in (b) and (d) 
shows bacltground estimated from the H scans. Shaded areas, 
which serve as guides to the eyes, show the spin-wave scatter- 
ing (yellow), and the extra scattering (gray), (e) and (f), and 
(g) and (h) show constant-energy scans at 30 meV around 
(1,0,5) and (3,0,1), respectively. The inset in (g) shows a H 
scan at L = 0. Black solid lines in (e) and (g) are guides to 
the eyes. Red squares in (f) indicate background measured 
around (1.2,0,5). The rise on the low-g side is due to the main 
beam. The red dotted line in (h) shows the spin- wave scat- 
tering with the above parameters. All data were measured 
with a counting time of 30 min or higher, but are normalized 
to counts/30min. 



scattering prefactor A{huj) at (1,0,5) and (3,0,1) can be 
explained by the transverse spin waves given the ordered 
moment along a. A single high point in Fig.^b) is spuri- 
ous due to neutrons scattered off the aluminum frame of 
the analyzer. It was not observed at any other antiferro- 
magnetic wave vectors or energies and hence has a negli- 
gible effect on Mj. At hu; = 30 meV [Figs.HJe) and 4(f)), 
the spin-wave scattering disappears or becomes negligi- 



bly small and the magnetic excitations are dominated by 
the extra scattering. Figures [3Je) and 4(g) show sharp 
peaks along H around (1,0,5) and (3,0,1), which are also 
observed at a different L as shown in the inset for L — 0. 
On the other hand, Figs. [Iff) and 4(h) show very broad 
scattering along L. As shown in Fig. IHh), this broad- 
ening along L cannot be accounted for by the spin-wave 
scattering. This extra scattering exhibits an anisotropic 
character with the broad L dependence and sharp H de- 
pendence, suggesting two-dimensional spin dynamics. 

This anisotropic scattering at high energy is possibly 
indicative of electron-hole excitations inside the Stoner 
continuum. In Fig. [S^a), we observe the rise of S(huj) 
starting at 15 meV and reaching a maximum around 
20 meV (shown by the gray shaded area), which de- 
viates from the l/huj dependence shown by the dotted 
line. According to a weak-coupling theory;^ an en- 
ergy gap, above which the collective spin-wave excita- 
tions are expected to merge into the Stoner electron- 
hole continuum, is given by = l.TGa/csTiv, where 
a — {vi + V2)/ {4:ViV2)^^'^', vi and V2 are Fermi veloci- 
ties of electrons and holes, respectively. For BaFe2As2, 
vi ~ V2 ~ 0.5 eV A^ and therefore ~ 20 meV. 
As a comparison, assuming the same Fermi velocities, 
As's are 26 and 33 meV for CaFe2As2 and SrFe2As2, 
respectively. The larger As's might explain why this 
anisotropic scattering has not yet been observed in 
CaFe2As2 and SrFe2As2»^'^ Similar broad inelastic- 
scattering lobes, which represent magnetic excitations 
inside the Stoner continuum, have been observed in 
V2-y03 (As = 1.4 meV) around incommensurate Bragg 
peaks below Tjv^ In BaFe2As2, we possibly observe 
the coexistence of the collective spin-wave excitations 
and electron-hole excitations near the edge of the Stoner 
continuum at 24 meV and only the single-particle spin 
dynamics above 30 meV. However, the presence of the 
spin-wave scattering inside the Stoner continuum and the 
anisotropic character appear to differ from V2-y03. Fur- 
ther experimental work on BaFe2As2, as well as on other 
iron arsenide compounds, is needed to further study these 
two-dimensional spin dynamics at higher energy. 



VI. CONCLUSION 

The magnetic excitations in BaFe2As2 have been mea- 
sured using inelastic neutron scattering. We observe low 
energy spin-wave excitations with a spin gap in the or- 
dered state. The smaller ratio of the out-of-plane to 
in-plane spin-wave velocities compared to CaFe2As2 and 
SrFe2As2 suggests that the magnetism in BaFe2As2 is 
more two-dimensional. At high energy, we observe the 
anisotropic scattering centered at the antiferromagnetic 
wave vectors. This scattering possibly indicates two- 
dimensional spin dynamics inside the Stoner electron- 
hole continuum. In the vicinity of T/v, the spin gap 
closes, and the spin fluctuations are observed up to 12 
meV (the highest energy of our measurements). Above 
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T/v, these spin fluctuations show the anisotropic charac- 
ter, attesting to the two-dimensionahty of this system. It 
is beheved that two-dimensionahty plays a significant role 
in high- Tc superconductivity. The fact that potassium- 
doped BaFe2As2 has the highest Tc so far among the 
nonoxide iron arsenide compounds hints at the con- 
nection between superconductivity and two-dimensional 
magnetism. It is interesting to see if two-dimensional 
spin fluctuations exist in the superconducting state. 
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